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We study the transitions between neighboring energy levels in a quasi-one-dimensional semicon- 
ductor quantum dot with two interacting electrons in it, when it is subject to a linearly time- 
dependent electric field. We analyze the applicability of simple two-level Landau-Zener model to 
describe the evolution of the probability amplitudes in this realistic system. We show that the 
Landau-Zener model works very well when it is viewed in the adibatic basis, but it is not as robust 
in the diabatic basis. 
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I. INTRODUCTION 



The Landau-Zener (LZ) model describes in a simplified way the ubiquitous situation of avoided crossings of energy 
levels in quantum mechanics This happens in numerous areas of physics such as quantum optics, atomic physics, 
nuclear physics, etc. In spite of its simplified nature, the LZ model often captures the essential features of avoided 
level crossings in realistic systems and the list of its applications continues to grow. 

Recently, we have proposed a generic quantum control method based on the navigation of the energy spectrum. 
D | The navigation of the spectrum is done varying a control parameter diabatically and adiabatically. The possibility 
^ . of traveling through a complex spectrum depends crucially on the nature of the energy level crossings. Our method 
requires that the system behave locally (at avoided crossings) like a LZ model, in the sense that complete diabatic 
and adiabatic transitions be possible. So, the first step in the application of this method to a realistic system must be 
a careful examination of the validity of this condition. Note that in a realistic system the interaction between levels 
is often intricate and the possibility that the LZ model worked has been discussed 0, Q . Recently, we successfully 
I ■ applied this control strategy to a quantum dot system and the isomerization of a LiCN molecule. In the present 
paper we analyse in detail the issue of the applicability of the LZ model to the avoided crossings of the two-electron 
quantum-dot system studied in Refs. [1, Q. 

Quantum dots are prime candidates to study the ideas and proposals of quantum control, given their flexible and 
tunable properties. In this paper, we continue the study of a quasi-one-dimensional double-dot system with two 
interacting electrons. Because of the one-dimensionality, this system is well suited to investigate new methods of 
quantum control, and at the same time it incorporates the important aspect of the interparticle interaction treated 
• exactly. We remark that the presence of interactions between particles is crucial in the new science of quantum 
[ information processing. 

\Q ■ The article is organized as follows: In order to make the work as self-contained as possible, in the next section we 
review the well-known Landau-Zener Model. In Section IIIII we describe our system and in the following section we 
, present the results concerning the applicability of the LZ model. We conclude with some final remarks. 
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O ■ II. LANDAU-ZENER MODEL 



The LZ model [l[ attempts to describe the universal situation of two levels interacting at an avoided crossing when 
a parameter A in the Hamiltonian is varied. The model consists of a two-level system described by a parameter- 
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dependent Hamiltonian, which expressed in the diabatic basis reads 
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(1) 



where 5 is a constant while e\ and £2 are linear functions of A: e\ = e + «i(A — A), £2 = £ + «2(A — A). The center 
of the avoided crossing is located at A and e [see Fig. [Q. The diabatic basis, |1) and |2), are parameter-independent 
eigenstates of the Hamiltonian Eq. (fTJ) with (5 = 0. 

The eigenenergies £'i(A) and i?2(A) of the Hamiltonian ((T|) are two hyperbolas (the adiabatic curves) as shown 
in Fig. [T] The eigenstates associated to those energies are the so-called adiabatic states, which we denote |</>i(A)) 
and |</> 2 (A)). The asymptotes to the energy hyperbolas are the diabatic straight lines £i(A) and £2(A). The shortest 
distance between the hyperbolas is 26. 
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FIG. 1: Schematic plot of an avoided crossing. The solid lines represent the adiabatic levels -Ei(A) and E2 (A). The asymptotes 
to the energy hyperbolas are the diabatic straight lines £i(A) and 62(A) (dashed lines). The center of the avoided crossing at A 
and e is inidicated with dotted lines. 



In his seminal paper, Zener considered the parameter A as a linear function of time and obtained the asymptotic 
probabilities of transitions between the diabatic states in this time-dependent problem. Assuming that the state |1) is 
the initial state (at t — > —00) and \{t) = /3t, and calling \ip(t)} the evolving wave function, the asymptotic probability 
to end up in the other diabatic state is 



00) = \{2\ip(t -» oo))| 2 = 1 - exp 
-2ttS 2 
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The two limiting cases in terms of (3 (the rate of change of the parameter A) are: 
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In this case the system follows the adiabatic curve going from the initial diabatic 
-r. The evolution takes place on the diabatic curve and the system remains in the 
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i) Slow transition: (3 <C 
state to the other one. 

ii) Rapid transition: [3 3> h \ a 
initial diabatic state. 

These limiting cases play a central role in our control method Q , as they give us quantitative criteria to choose either 
the diabatic or adiabatic paths in traversing an avoided crossing 

More recently, the complete time dependence of the occupation probabilities in both the diabatic (|(l|^;(f))| 2 , 
\{2\ip(t))\ 2 ) and adiabatic basis (| (0i(A) |-0(^)> | 2 , | <0 2 (A) |V(*)> | 2 ) sets has been obtained @,0]. 



III. THE SYSTEM: QUASI-ONE-DIMENSIONAL DOBLE QUANTUM DOT WITH TWO 

INTERACTING ELECTRONS 



Let us consider a quasi-one-dimensional double quantum dot with two interacting electrons in the presence of a uni- 
form longitudinal electric field. This system is interesting for two reasons. First, this type of system is experimentally 
realizable nowadays, and second, the nonperturbative interparticle interaction is taken into account. 

We have chosen a semiconductor system with realistic dimensions: it is an elongated quantum dot 100 nm long and 
50 A wide. Due to the small thickness of the structure, the energies of the transverse modes are widely spaced and it 
is enough to consider only the lowest transverse state. Therefore, an effective Hamiltonian that depends only on the 
longitudinal coordinate z describes the dynamics of the system [8[ 

ft 2 d 2 d 2 

H ^-^d4 + d4 ) + V{zi) + V{Z2) 

+V c (\zi - z 2 \) - e( Zl + z 2 )E{t) , (3) 

where m is the electron effective mass in the semiconductor material, Vc is the Coulomb interaction between the 
electrons, V is the confining potential, and E(t) is a time-dependent external electric field. There is no restriction 
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FIG. 2: Confining double well potential in the longitudinal direction of the coupled quantum dot structure. The external 
electric field is E — (solid lines) and E — 12kV/cm (dashed lines). 



regarding the choice of the confining potential in the z-direction, but we have selected the double-well configuration 
shown in Fig. [51 Double-well potentials are interesting due to the interplay between tunneling and localization. 
Moreover, in our two-electron system these important phenomena can be related to the Coulomb interaction. In all of 
the time evolutions that we will analyze, we assume that the wave function is initially a singlet (antisymmetric spin 
wave function). Since the Hamiltonian is spin independent, the spin wave function remains a singlet and the orbital 
part of the wave function is symmetric at all times. 

We have used the time-dependent electric field E(t) as the control parameter The first step is to understand 

the behavior of our system when the electric field is taken as a constant. So, we have computed numerically the 
eigenergies and eigenfunctions of the system as a function of a constant electric field. We numerically diagonalize the 
Hamiltonian of Eq. §3§ expanded in the basis set of Slater determinants constructed with the 12 bound single-particle 
states of the double-well potential. The two-particle basis set has then 12*(12+l)/2=78 states. 

The energy spectrum as a function of the external electric field is shown in Fig. IIII1 We can see that the spectrum 
is composed by fairly straight lines which never cross each other, resulting in frequent avoided crossings. For the low- 
lying states included in our spectrum, and far from the avoided crossings, the adiabatic states have clear localization 
properties connected with their slope (see Fig. 2 in Ref. [H): 

i) in the eigenstates corresponding to negative slope both electrons are in the left well, 

ii) the states associated to the positive slope are localized on the right well, and 

iii) the states with neutral slope have one electron in each well. 

At avoided crossings, the eigenstates mix their localization characteristics reaching the maximal degree of mixing at 
the center of the avoided crossing. 



IV. ANALYSIS OF THE APPLICABILITY OF THE LANDAU-ZENER MODEL 



In ourprevious works we introduced a method of quantum control via traveling in the energy spectrum of a quantum 
system [H, [||. The building blocks of this method are, on the one hand, the adiabatic evolutions far from avoided 
crossings, and on the other, the slow and fast evolutions employed at avoided crossings in order to shift in a controlled 
way from one adiabatic path to a neighboring one. If the system behaves locally like a LZ this possibility will be 
guaranteed. For this reason, in this section we will analyze the range of validity or applicability of the Landau-Zener 
model to describe the transitions at the avoided crossings of our system. 

We begin by studying the avoided crossing labelled "A" in Fig. IIIIl between the ground state and the first ex- 
cited state near the value of the electric field E = 5kV/cm. Initially the system is in the ground state with no 
electric field (state labelled "a" in Fig. IIII[) and we study the probability to remain in the ground state when the 
electric field is increased linearly with time at different velocities. This corresponds to the adiabatic probability 
\(4>i(E)\ij){t))\ 2 introduced at the end of Section ITT] These probabilities are shown in Fig. 2J We remark that in 
order to compare the results for different velocities, we plot in Fig. 2] the probabilities as functions of the electric 
field rather than as functions of time. In Fig. [1] we present the adiabatic probability for the following velocities, 
E = 0.07,0.27,0.53, 1.07, 4.27kV/cmps. 

We now compute the adiabatic probabilities in the LZ model. The first step is to fit the parameters 6, the location 
of the avoided crossing, oti, and a-i [9( of the two- level Hamiltonian of Eq. JT]) to the avoided crossing under study. As 
initial state we take the adiabatic LZ state for the value of A that corresponds to state "a" in Fig. IIIIl We compute the 
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FIG. 3: The energy spectrum of the two interacting electrons confined in a quasi-one-dimensional double-well semiconductor 
nanostructure as a function of an external uniform electric field. It is plotted the first 31 energy levels. See text for details. 



adiabatic probabilities in the LZ model for the previous set of rates of change of the electric field. These results are 
plotted with solid lines in Fig. Q] We can see that our system is well described by the LZ model for the whole range 
of velocities considered. Moreover, we see that the asymptotic probability obtained by Zener (horizontal dashed lines 
in Fig. 0| gives accurate results far from the avoided crossing. Note that far from the avoided crossing the diabatic 
and adiabatic states in the LZ model are essentially the same. 

We have done similar analyses for other avoided crossings of the energy spectrum. Namely, we start with the states 
labelled "b", "c", and "d" in Fig. IIII1 and we study the transition probabilities in the adjacent avoided crossings "B", 
"C", and "D", respectively. In Fig. [5] we show the results, which verify the previous conclusion, in the sense that in 
the adiabatic basis the two-level LZ model fits very well the exact results. It is worth noting here that in our previous 
work of Ref. 0, [f| we travel in the spectrum (that is, we attempt to go from a given adiabatic state to another one). 
In this sense, the above given results are the most relevant ones to judge the applicability of our control method. 

Since the LZ model is defined on the basis of diabatic states, it is perhaps more natural to perform the former 
analysis on that basis set. However, the question arises of what the diabatic states are in our realistic system. Indeed, 
in a multilevel system like ours the two states involved in the avoided crossing become mixed with other states and 
therefore acquire a dependence on the control parameter (which is not allowed for in the usual LZ model). It is thus 
an interesting question to ask whether it is possible to find a "fixed" basis set which could play the role of the diabatic 
basis in the LZ model. For example, we now calculate the probability \(4>i(Eo)\i{j(t))\ 2 , where |<^i(i?o)} is the the 
initial state in the dynamic passage of an avoided crossing. That is, we are considering \4>i(Eq)) as being one of the 
diabatic basis states. We now do this for the lowest crossing taking Eq = 0, and compare with the results of using 
the two-level LZ model in Fig. [6] One can clearly see that the agreement between the two calculations is not very 
good. This can be understood with the help of the inset of Fig. [HI which shows the overlaps \{<p\(E = 0)\(j>i(E))\ 2 and 
\{4>i(E — O)|02(£')}| 2 as functions of the electric field E. It is clear from the inset, especially from \{4>i{E — 0)\4>2(E))\ 2 
(dashed line), that the hypothesis of a parameter- independent diabatic state is not satisfied (that the overlap is not 
equal to one far from the avoided crossing), and therefore the LZ model tends to fail. However, in other avoided 
crossings we have observed that it is possible to find good diabatic states (which are fairly parameter-independent 
around the avoided crossing). For example, we repeated the previous analysis for the avoided crossing labelled 
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FIG. 4: Adiabatic transition probabilities \{<f>i(E)\rp(t)}\ 2 at the lowest avoided crossing (labelled "A" in Fig. IHIjl for various 
rates of change of the control parameter. The velocities are: E = 0.07(A), 0.27(D), 0.53(x), 1.07(o), 4.27 (*) kV/cm ps. The 
initial state in the exact evolution is the one labelled "a" in Fig. Mil The solid lines give the adiabatic probabilities in the LZ 
model and the dotted lines are the asymptotic LZ probabilities given in Eq. |2jl. 




Electric field (kv/cm) 



8 10 12 14 

Electric field (kv/cm) 



7.5 8 8.5 9 9.5 1 

Electric field (kv/cm) 



FIG. 5: Adiabatic transition probabilities at various avoided crossings (labeled "B", "C" and "D" in Fig. IIII[) . The initial 
states in these evolutions are labelled in Fig. IHII as "b" [panel (a)] , "c" [panel (b)], and "d" [panel (c)] . The velocities are: (a) 
E = 0.0015(A), 0.0077(D), 0.015(x), 0.07(c), 4.27 (*)kV/cm ps, (b) E = 0.07 (A), 0.53 (D), 1.27 (x), 4.27 (o), 40.0 (*) kV/cm ps, 
and (c) E = 0.003 (A), 0.07 (D), 0.27 (x), 0.53 (o), 4.27 (*) kV/cm ps, 



"B" in Fig. IIIII choosing \<fiiQ(E — 1.4kV/cm)) as one of the diabatic states. In Fig. we plot the probability 
\((pxe(E = 1.4 kV/cm) | 2 , which shows a better agreement than the one in Fig. [6j We remark that, as can 
be seen in the inset of Fig. [7J the state \4>\q{E = 1.4kV/cm)) is a good choice of diabatic state, since the overlap 
\(<f>ie(E = 1.4kV/cm)|</>i7(.E))| 2 is close to one at the right of the crossing and close to zero to the left (see dashed 
line). The behavior seen in the inset is exactly what one obtains in the LZ model for the overlaps between the diabatic 
and adiabatic bases. 



V. FINAL REMARKS 



We have studied the applicability of the LZ model in a realistic system: a quasi-one-dimensional double quantum 
dot with two interacting electrons. We showed that the LZ model works very well when it is viewed in the adiabatic 
basis. This result is the cornerstone for the quantum controlability using the method of control introduced in 

However, when seen in the diabatic basis the results are not so robust as in the case of the adiabatic basis. This 



6 




FIG. 6: Diabatic transition probabilities \{(f>i(Eo = 0)\ip(t))\ 2 (dotted lines) at the lowest avoided crossing, labelled "A" in Fig. 
Mil The rates of change of the control parameter, the electric field E, are the same as Fig. [4] In solid lines, the corresponding 
probabilities in the LZ model. Inset: overlaps \(</>i(E = O)|0i(i?))| 2 (solid line) and \{<f>i(E = 0)\()>2(E)}\ 2 (dashed line). 
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FIG. 7: Diabatic transition probabilities \(<j>ie{E = 1.4kV/cm)[*0(t)}[ 2 (dotted lines) at the avoided crossing, labelled "B" 
in Fig. IIIII The rates of change of the control parameter, the electric field E, are the same as in Fig. [5ja). In solid lines, 
the corresponding probabilities in the LZ model. Inset: overlaps \{(f>ie(E = 1.4kV/cm)|</!>i6(-E))| 2 (solid line) and \(<J>i6(E = 
1.4kV/cm)|(£i7(£))| 2 (dashed line). 



is due to the fact that, for multilevel systems, a proper diabatic basis does not exist. Rather, the pair of interacting 
levels at an avoided crossing become mixed with other states and acquire a dependence with the control parameter 
even far from the avoided crossings. 
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